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Abstract 

A Hadamard matrix of order re has no proper Hadamard subma- 
trices of order m > re/2. We generalise this result to maximal deter- 
minant submatrices of Hadamard matrices, and show that an interval 
of length ~ re/2 is excluded from the allowable orders. We make a 
conjecture regarding a lower bound for sums of squares of minors of 
maximal determinant matrices, and give evidence to support it. We 
give tables of the values taken by the minors of all maximal deter- 
minant matrices of orders < 21 and make some observations on the 
data. Finally, we describe the algorithms that were used to compute 
the tables. 



1 Introduction 



A {+1, — l}-matrix (abbreviated "{±l}-matrix" below) is a matrix whose 
elements are +1 or —1. The Hadamard maximal determinant problem, posed 



by Hadamard [T5], is to find the maximal determinant D(n) of an n x n 
{±l}-matrix of given order n. A matrix that attains the maximum is a 
maximal determinant matrix (abbreviated maxdet matrix). Such matrices, of 
which Hadamard matrices are a special case, are of interest in combinatorics 
and have applications in statistical design (20], coding theory and signal 
processing (TTJ [321 EE]- m design theory a maxdet matrix is also known as a 
saturated D-optimal design. 

Finding maxdet matrices is in general difficult, with the difficulty depend- 
ing on the congruence class of n (mod 4). Most is known for n = (mod 4) 
(the Hadamard orders), and least for n = 3 (mod 4). Maxdet matrices are 
known for orders 1 through 21 inclusive; for n = 22 and n = 23 there are 
conjectures (see [52]) but as yet no proofs of maximality. Constructions exist 
for various infinite families and ad- hoc examples, see [33J. 

Consider the k x k submatrices Sk(A) of an n x n {±l}-matrix A, where 
< k < n. For M £ Sk{A), we say that det(M) is a minor of order k of A. 
Usually only the absolute value of det(M) is of interest. 

One method for finding maxdet or large-determinant matrices involves 
choosing a large-determinant matrix of a slightly smaller (or larger) order 
than the desired order, possibly perturbing it by a low-rank modification, and 
adding (or removing) a small number of suitably chosen rows and columns. 
For example, Solomon [32, 39J found a (conjectured) maxdet matrix of order 
33 and determinant 441 x 2 74 in this way by starting with an appropriate 
(large excess) Hadamard matrix of order 32, and in [SJ EH EE] the method 
was used to obtain lower bounds on D(n). 

This motivates our interest in the minors of maxdet matrices, and in par- 
ticular the question: What is the largest order of a maxdet matrix contained 
as a proper submatrix of a given maxdet matrix? In this paper we answer 
this question for the maxdet submatrices of maxdet matrices of orders n < 21 
by computing all minors of maxdet matrices of these orders. Our work ex- 
tends that of earlier researchers who have considered minors of Hadamard 
and maxdet matrices with other applications in mind, such as the problem 
of growth in Gaussian elimination [TT1 E3J ESI EE1 133] . 

Schmidt [541 pg. 441] says "The nature of the construction ■ ■ ■ is in line 
with the computer assisted observation that binary matrices with maximal de- 
terminants may not contain large order submatrices with large determinants". 
Whether this is true in general depends on the precise meaning of "large". 
Certainly there are exceptions. For example, the maxdet {±l}-matrix of 
order 17 contains a maxdet (Hadamard) submatrix of order 16. 
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In Theorem [T] of §2] we prove that a Hadamard matrix of order n can 
not have a proper Hadamard submatrix of order m > n/2. We then con- 
sider maxdet submatrices of Hadamard matrices, and show, assuming the 
Hadamard conjecture, that a Hadamard matrix of order n can not have a 
maxdet submatrix M of order m > n/2 + 5 logn unless m > n — 2 (see The- 
orem [2]). In other words, m can not lie in the interval (n/2 + 5 logn, n — 2). 
Without the assumption of the Hadamard conjecture, we can still exclude 
an interval (n/2 + o(n),n — o(n)) of length ~ n/2, see Theorem [3] and Re- 
mark [2j These results partially confirm the remark of Schmidt quoted above. 
However, they apply only to sub-matrices of Hadamard matrices. We can 
not, except in small cases that are amenable to explicit computation, exclude 
the possibility that a maxdet matrix of order n ^ (mod 4) has a maxdet 
submatrix of any given order m < n. 

In §|3] we define two sequences related to the sets of minors of maxdet 
matrices, and give the first 21 terms of each sequence. 

In §4] we describe the minors that occur in maxdet matrices of orders 1 
through 21, and make some observations on the patterns that occur. We 
mention a result (Proposition [3]), on small minors of Hadamard matrices, 
which was suggested by the data before a proof was found. 

Motivated by Turan's result [U] that the expected value of det(A) 2 is 
m\ for random {±1} matrices of order m, we consider the mean value of 
det(M) 2 over all m x m submatrices of maxdet matrices of order n > m, and 
conjecture that it is bounded below by m\ (see Conjectured]). The conjecture 
is consistent with our computations for 1 < m < n < 21. 

Finally, in §|5] we describe the algorithms that we used to compute mi- 
nors of square {±l}-matrices, as well as some related algorithms that were 
considered but rejected for various reasons. 

Hadamard equivalence and HT-equivalencefj] 

Two {±l}-matrices are Hadamard equivalent if one can be obtained from the 
other by negating rows or columns, and/or by interchanging rows or columns. 
Clearly the answer to the question posed above is the same for all matrices in 
a Hadamard equivalence class, and also for any matrix A and its transpose 
(dual) A T . Hence, it is useful to define the notion of HT- equivalence by saying 

lr There seems to be no widely- accepted name for this concept. In [5] HT-equivalence 
is called "extended Hadamard equivalence" . Wanless [42. calls an HT-equivalence class a 
"resemblance class" . 



3 



that two matrices A and B are HT- equivalent if A is Hadamard-equivalent 
to B or A is Hadamard-equivalent to B T . For example, there are 5 distinct 
Hadamard equivalence classes of Hadamard matrices of order 16, but two of 
these classes contain matrices that are dual to those in the other class, so 
there are only 4 distinct HT-equivalence classes. 

2 Excluded minors of Hadamard matrices 

In this section we consider the possible orders of submatrices M of a Hadamard 
matrix H, satisfying the condition that M is Hadamard (see Theorem [1]) or, 
more generally, that M is a maxdet matrix (see Theorems I2H3]) . 

Recently Szollosi [101 Proposition 5.5] established an elegant correspon- 
dence between the minors of order m and of order n — m of a Hadamard ma- 
trix of order n. His result applies to complex Hadamard matrices, of which 
{±1} Hadamard matrices are a special case. More precisely, if d + m = n, 
< d < n, then for each minor of order m and determinant /x there corre- 
sponds a minor of order d and determinant ±n n / 2 ~ d /i. Only a few special 
cases (for small d or m) were known before Szollosi's result (see for ex- 
ample [TTJ |25l E51 EZ]), although we note that Szollosi's Lemma 5.7 follows 
easily from Jacobi's determinant identity [18]. We use the following corollary 
of Szollosi's theorem. 

Corollary 1. Suppose that a Hadamard matrix H of order n has an m x m 
submatrix M , where n/2 < m < n. Then 

|det(M)| < n m - n/2 D{n - m) . 

Proof. By Szollosi's theorem, | det(M)| = n n / 2 ~ d \ det(M')|, where d = n — m 
and M' is some d x d submatrix of H. Since n/2 — d = m — n/2 and 
| det(M')| < D(d) by the definition of D, the corollary follows. □ 

The following lower bound on D(m) is given in [8j Corollary 3]. 

Proposition 1. Assume the Hadamard conjectur^. Then, for m > A, we 
have D(m) > 4m m / 2 ~ 1 . 

2 It is sufficient to assume that Hadamard matrices of order 4k exist for all positive 
integers k < (m + 2)/4. This is known to be true for Ak < 668, see [?T1 136] . 
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Suppose that a Hadamard matrix H of order n has a maxdet submatrix M 
of order m. Corollary [T] gives an upper bound on | det(M)|, and Proposition [TJ 
gives a lower bound. Theorems [TH2] show that these bounds are incompatible 
for certain values of m. Theorem [1] considers the case that M is a Hadamard 
matrix, and Theorem [2] considers the more general case that M is a maxdet 
matrix. We have stated Theorem [2] under the assumption of the Hadamard 
conjecture, but a weaker result is provable without this assumption - see 
Theorem |3j 

Theorem 1. Let H be a Hadamard matrix of order n having a Hadamard 
submatrix M of order m < n. Then m < n/2. 

Proof. The theorem is trivial if m < so assume that n > m > n/2. 
Let d = n — m. Since M exists and | det(M)| = m m l 2 , Corollary [1] and 
Hadamard's bound for D(d) give 

m m/2 < n m-n/2 d d/2 _ ^ 

Squaring both sides of (JT]), we have (m/n) m < (d/n) d . Taking n-th roots 
and defining x := m/n G (0, 1), we see that 

x x < (l-xY~ x . 

This inequality is equivalent to f(x) < 0, where / : [0, 1] — > M is defined by 

J xlnx — (1 — x) ln(l — x) if x £ (0,1), 
f(x) — < 

I otherwise. 

It is easy to verify that /(1/2) = /(l) = 0, f'{x) = —2 — In x — ln(l — x), and 

/ "w = ^)< 0ta ( 1 / 2 - 1 »- 

Thus, /(x) > in (1/2,1), so we must have x < 1/2 or x = 1. The case 
x = 1 is ruled out because it implies that m = n, contrary to the assumption 
that m < n. Thus x < 1/2, which implies that m < n/2. □ 

Theorem 2. Assume the Hadamard conjecture. Let H be a Hadamard 
matrix of order n having a maxdet submatrix M of order m < n. Then 
Tn < (n/2 + 5 Inn) or m > n — 2. 
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Proof. The result is trivial if m < n/2, and n/2 + 51nn>n — 3 for n < 28, 
so we can assume that m > n/2 > 14. By Proposition [T] and Corollary [TJ we 
have 

m m/2 (±\< n m~n/2 d d/2^ (2) 

\m J 

where d = n — m. We prefer to use the slightly weaker inequality 

m m/2 f±\ < n m-n/2 d d/2_ (3) 

Taking logarithms, and defining x := m/n and / as in the proof of Theorem [H 
we obtain 

/<«) < (4) 
n 

The right side of this inequality is positive (since n > 4) and independent 
of x (this is why we used (J3J) instead of (j2J)). We showed in the proof of 
Theorem H] that /(1/2) = /(l) = and f (x) < on (1/2,1). Let 



X -i j 



(1 + v 7 ! - 4/e 2 )/2 w 0.84 



be the (unique) point in (1/2, 1) at which f'(x) vanishes. Thus f(x) attains 
its maximum value at x — x max . Since 21n(n/4)/n < 0.14 < f(x max ) ~ 0.15 
for n > 28, the inequality (j3]) is not satisfied for all x G (1/2, 1), and there 
is a unique interval (xo,Xi) C (1/2,1) on which f(x) > ln(n/4)/n, with 
Xmax ^ (^Oj^i)- (Here x and x\ depend on n but not on m.) It follows that 
there can not exist a maxdet submatrix of order m with xq < m/n < X\. 
To locate X\ we consider the case d = n — m = 3. The inequality ([3]) gives 

n-3\ n ~ 3 27 

<77-' ( 5 ) 



n / 16n 

but the left side of this inequality is bounded away from zero as n — > oo, 
whereas the right side tends to zero. Thus, the inequality can not hold for 
large n. In fact, a computation shows that dSJ) can only hold for n < 29. 
Thus, for n > 29 an interval (xo,xi) as above exists, with xi > 1 — 3/n, so 
m = n — 3 is not a possible order of a maxdet submatrix M. 

We now show that nxo < n/2 + 51nn. Define 5 := m — v. Thus m = u+5, 
d = v — 5, and squaring the inequality (j3j) gives 

(^ + 5)" +5 0) < (2z/) 25 (z/-5)^, 
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or equivalently 

'iy + SY /i/\2 / 4z/ 2 



Write z := 8/v = 2x — 1. We can assume that 



<U kF« ■ ( 6 ) 



< 2 < z maa; = 2ac max - 1 = y/l - 4 /e 2 ~ 0.68 . 
Taking logarithms in flS]) gives 



-In I < <^(ln4-ln(l -2 2 )) + 21n(z//2). (7) 

Collecting the terms involving 5 gives 

5(2-ln4-£(») < 21n(z//2), 

where 

e(z)=2- 



i hl (I±£)_ ln(1 _^ E 

v 7 fc=l 



z 2k 



z \l-zj ' ^k(2k + l) 

is monotonic increasing on [0, z m3X ], so 

e{z) < e(z max ) < 0.1803, 

and thus 2 - In 4 - e(z) > 0.4334 for z G (0, z max ). It follows that 

. 21n(i//2) 

5 < y —L_L < 5 In n. 

0.4334 



□ 



Remark 1. We do not expect the values m = n- lorm = n- 2to occur 
for n > 4. They have to be included as possibilities simply because the 
lower bound on D(ri) given by Proposition [1] is too weak to exclude them. 
It is possible to have m > n/2, for example a maxdet submatrix of order 
m = n/2 + 1 occurs in Hadamard matrices of orders n = 4 and n = 12. 

We can prove a result similar to, but weaker than, Theorem [2] without 
assuming the Hadamard conjecture. Let the prime gap function X(n) be the 
maximum gap between consecutive primes (pi,Pi+i) with pi < n. Then we 
have: 
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Theorem 3. Let H be a Hadamard matrix of order n > 4 having a maxdet 
submatrix M of order m < n, and let X(n) be the prime gap function. 
Then there exist positive constants C\, C2 such that m < n/2 + ciA(n)lnn 
or m > n — C2\{n). 

Sketch of proof. The proof is similar to that of Theorem |2J but uses Theo- 
rem 1 and Corollary 1 of j8] in place of Proposition [TJ Thus ([3]) is replaced 
by 

/ , n A(n/2)/2 
m m/2 M < n m-n/2 dd /2 ; 

\ne J 

and gD by 

ft n ^ w / Q x ln(we/4) 
/(x) < A(n/2) 

The remainder of the proof follows that of Theorem [2j except that some 
of the explicit constants have to be replaced by 0(A(n/2)) terms, and we 
have to assume that n is sufficiently large, say n > n . At the end, we can 
increase c\ or C2 if necessary to ensure that n/2 + ciA(n) Inn > n — C2\{n) 
for 4 < n < hq. □ 

Remark 2. By a result of Baker, Harman and Pintz [2], \{n) = 0(n 21 / 40 ), 
so the excluded interval in Theorem [3] has length ~ n/2. 



3 Sequences related to minors 

In this section we define two sequences related to the minors of maxdet 
matrices, and give the first 21 terms in each sequence |31j . 

For the convenience of the reader, Tables 1-2 give some data taken from 
Orrick and Solomon [32], where references to the original sources may be 
found. Table 1 gives the spectrum of possible (absolute values of) deter- 
minants of {±l}-matrices of order n < 11, normalised by the usual fac- 
tor 2 n_1 . In this and other tables, the notation "a.. b" is a shorthand for 
ll {x E N : a < x < 6}". Table 2 gives A(n) := D(n)/2 n ~ 1 for n < 21. 

We are interested in when the full spectrum of possible minor values 
occurs in the minors of maxdet matrices of given order n. In the following 
definitions, N denotes the positive integers. 
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n 


Spectrum { | det( A) \/2 n 1 | 


1 


I 1 } 


2 


{0,1} 


Q 
O 


fn ii 


4 


{0-2} 


c 



{0..3} 


6 


{0..5} 


7 

8 


{0..9} 

{0..18,20,24,32} 


9 


{0.-40,42,44,45,48,56} 


10 


{0.102,104, 105, 108, 110, 112, 


116,117, 120,125,128,144} 


11 


{0..268, 270..276, 278..280, 282..286, 


288, 291, 294..297, 304, 312, 315, 320} 



Table 1: Spectrum of {±l}-matrices of order n < 11, 
from Orrick and Solomon [32J; for n = 13 see [7|. 

Definition 1. T/ie full-spectrum threshold o/ an n x n {±1} matrix A is 
the maximum rrif < n such that the full spectrum of possible values occurs 
for the minors of order rrif of A. 

Definition 2. The full-spectrum threshold rrif : N — > N is the maximum of 
the full-spectrum threshold of A over all maxdet matrices A of order n. 

We write rrif {A) or rrif(n) to denote the full-spectrum thresholds of Def- 
initions [T] or [2] respectively; which is meant should be clear from the context. 
The values of rrif(n) for 1 < n < 21 are given in Table 3. We note that 
the full-spectrum threshold rrif {A) does depend on the HT-equivalence class 
of A. For example, the four HT-equivalence classes for order 16 give four 
different values rrif G {5, 6, 7, 8}, see Tables 23-26. 

For the reasons mentioned in $11 we are also interested in the largest 
order of a maxdet matrix contained as a proper submatrix of a given maxdet 
matrix. We make some definitions analogous to Definitions [[HS 

Definition 3. The complementary depth of an n x n {±1} matrix A is the 
maximum < n such that a maxdet matrix of order occurs as a proper 
submatrix of A, or if n = 1. The depth of A is d(A) := n — m d (A). 



9 



n 


A(n) 


n 


A(n) 


n 


A(n) 


n 


A(n) 






1 


1 


2 


1 


3 


1 


4 


2 


5 


3 


6 


5 


7 


9 


8 


4 x 2 3 


9 


7x 2 3 


10 


18 x 2 3 


11 


40 x 2 3 


12 


6 x 3 5 


13 


15 x 3 5 


14 


39 x 3 5 


15 


105 x 3 5 


16 


8 x 4 7 


17 


20 x 4 7 


18 


68 x 4 7 


19 


833 x 4 6 


20 


10 x 5 9 


21 


29 x 5 9 











Table 2: Maximal determinants A(n) = D(n)/2 n 1 , n < 21, 
from Orrick and Solomon |32|. 



Definition 4. T7ie complementary depth : N — > Z is the maximum of 
the complementary depth of A over all maxdet matrices A of order n. The 
depth d : N — > Z is defined by d(n) := n — m^n). 

We write d(A) or d(n) for the depths of Definitions [3] or [5] respectively; 
similarly for md(A) and m<i(n). Clearly d(A) depends on the HT-equivalence 
class of A- for example, see Tables 30-32 for the three HT-equivalence classes 
of order 18 with depths 7, 7 and 10. From Definition HI d(n) is the minimum 
of d(A) over all maxdet matrices A of order n, so d(18) = 7. Computed 
values of d(n), md(n) and rrif(n) for 1 < n < 21 are given in Table 3. It is 
clear from the definitions that mf(n) < m^n) for n > 1. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


d 


1 


1 


1 


1 


1 


1 


1 


4 


1 


2 


3 


5 


m d 





1 


2 


3 


4 


5 


6 


4 


8 


8 


8 


7 


m f 


1 


1 


2 


2 


3 


4 


6 


4 


6 


6 


7 


6 



n 


13 


14 


15 16 


17 


18 


19 


20 


21 


d 


6 


7 


8 8 


1 


7 


10 


10 


10 


m d 


7 


7 


7 8 


16 


11 


9 


10 


11 


m f 


7 


7 


7 8 


8 


8 


9 


8 


10 



Table 3: Depths d(n) of largest maxdet proper 
submatrices, complementary depth m^(n) = n — d(n), 
and full-spectrum threshold rrif(n); see Definitions [2HH 

If n = (mod 8) then d(n) = n/2 in the range of our computations. If 
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n = 4 (mod 8) then both d(n) = n/2 (for n = 20) and = n/2 — 1 (for 
n = 4, 12) are possible. We see from Table 3 that the computed values all 
satisfy d(n) < (n + l)/2. It is interesting that the value d — 1 occurs for 
n < 7, n = 9 and n — 17 (contrary to the remark of Schmidt quoted in §T]). 

4 Further results and observations on minors 

In Tables 6-43, which may be found in the Appendix at the end of the 
paper, we give computational results on the minors of maxdet matrices of 
order n < 21. Here we make some empirical observations on the results and 
mention a fact (Proposition [3]) that was suggested by them. 

Let k = \n/4\. Factors of the form k 2k , /c 2fc_1 , . . . , k 2 , k are present as 
we descend through the minors of maxdet matrices of even order n. For 
Hadamard matrices this is an easy consequence of the Hadamard bound 
n n / 2 and Szollosi's theorem, but for n = 2 (mod 4) we do not have a sim- 
ple explanation. Factors of the form k 2k ~ x , k 2k ~ 2 , . . . , k are present in the 
minors of maxdet matrices with order n = 1 (mod 4), and factors of the 
form k 2k ~ 2 , k 2k ~ 3 , . . . , k are present if n = 3 (mod 4). It is an open ques- 
tion whether this behaviour persists for n > 21. The observed divisibility 
properties are related to the structure of the Gram matrices A T A of maxdet 
matrices A, but in general this structure is unknown. For a summary of what 
is currently known, see [32J. 

The presence of high powers of k = [n/4:\ in the minors of order m, 
and high powers of a possibly different integer k' = [m/4\ in the maximal 
determinants of order m, gives one explanation of why certain minors can 
not meet the maximal determinant for that orderly For example, Table 19 
shows that a Hadamard matrix of order n = 12 has minors of order 11 with 
scaled value 3 5 , but A(ll) = 5 x 2 6 , which contains a high power of 2, not 
of 3. 

Proposition [2] is from [91 Theorem 1] . The upper bound is sharp because 
it is attained for Hadamard matrices. For the case that A is a Hadamard 
matrix, the result is due to de Launey and Levin [28, Proposition 2]. 



3 Of course, this begs the question of why the maximal determinants are divisible by a 
high power of k' - we do not have a convincing explanation for this unless the order m is 
such that the Hadamard, Barba [3], or Ehlich-Wojtas [T3l [44] bound is achieved, in which 
case it follows from the form of the relevant bound, see Osborn [33] [pp. 98-99]. 
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Proposition 2. Let A be a square {±1} matrix of order n > m > 1. Then 
the mean value E(det(M) 2 ) ofdet(M) 2 , taken over all m x m submatrices 
M of A, satisfies 



Moreover, equality holds in ^ iff A is a Hadamard matrix. 

Remark 3. For random {±1} matrices M of order m, the expected value 
of det(M) 2 is ml, by a result of Turan [4T]. Proposition [2] implies that 
£'(det(M) 2 ) > ml for the order- m submatrices M of a Hadamard matrix, 
with strict inequality if m > 1. 

As a check on the correctness of our programs, we computed the mean 
value of det(M) 2 for submatrices of order m of maxdet matrices of order 
n < 21, and 2 < m < n. The results agreed with the predictions of Proposi- 
tion [2j The following conjecture is consistent with our computations. 

Conjecture 1. Let A be a maxdet matrix. Then the mean value _E(det(M) 2 ) 
o/det(M) 2 taken over all mx m submatrices M of A satisfies the inequality 



Moreover, the inequality (Q is strict for m > 1 . 

Table 4 gives some data for orders 13 < n < 15 to support Conjecture [TJ 
In the table, Ri,(m,n) is the ratio of _E(det(M) 2 ), for submatrices M of 
order m of a maxdet matrix of order n, to the conjectured lower bound ml. 
Similarly, Rn(m,n) is the ratio of _E(det(M) 2 ) to the upper bound We 
see that R^im^) > 1 for 2 < m < n (as conjectured), and the lower bound 
is reasonably good for m < 5, but deteriorates for larger m. The upper 
bound is within a factor of three of _E((det M) 2 ) for all m. A similar pattern 
occurs for all orders n < 21, except that Rn{m,n) = 1 if n is a Hadamard 
order, in accordance with Proposition [2J 

The frequencies of occurrence of small singular submatrices of Hadamard 
matrices are given in the following Proposition (91 Corollary 4], which was 
suggested by the computational results before we found a proof. The case 
m = 2 is implicit in a paper of Little and Thuente [291 Pg- 254] . 





£(det(M) 2 ) > ml 



(9) 
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m 


R L (m, 13) R 


H (m, 13) 


R L (m, 14) R 


H (m,14) 


R L (m, 15) R 


H (m, 15) 


2 


1.077 


0.994 


1.067 


0.991 


1.059 


0.988 


3 


1.259 


0.983 


1.222 


0.973 


1.195 


0.966 


4 


1.611 


0.968 


1.516 


0.948 


1.445 


0.935 


5 


2.283 


0.949 


2.054 


0.917 


1.890 


0.897 


6 


3.625 


0.928 


3.072 


0.882 


2.694 


0.852 


7 


6.560 


0.904 


c 1 on 

5.139 


A O A O 

0.843 


4.233 


0.804 


8 


13.81 


0.879 


9.773 


0.802 


7.427 


0.752 


9 


34.80 


0.851 


21.58 


0.759 


14.79 


0.699 


10 


109.4 


0.823 


56.93 


0.715 


34.14 


0.645 


11 


457.4 


0.795 


187.0 


0.671 


94.00 


0.592 


12 


2864 


0.765 


815.6 


0.627 


321.7 


0.540 


13 


35796 


0.736 


5318 


0.584 


1461 


0.491 


14 






69137 


0.542 


9902 


0.444 


15 










133638 


0.399 



Table 4: Ratio of E(det(M) 2 ) to lower and upper bounds, 13 < n < 15. 
For definitions of R^m^n) and Ru(m,n), see text. 

Proposition 3. Let H be a Hadamard matrix of order n, and let Z(m, H) 
be the number of minors of order m of H that vanish. Then 

Z(2,H) = n 2 {n - l)(n - 2)/8, and (10) 

Z(3,H) = n 2 (n-l)(n-2)(n-4)(5n-4)/288. (11) 

Remark 4. There is no analogue of Proposition [3] for minors of order m > 3, 
because the value of Z(m, H) can depend on the HT-equivalence class of H, 
so is not given by a polynomial in m unless m < 3. For example, the 
four HT-equivalence classes of Hadamard matrices of order 16 have 1717520, 
1712912, 1710608, and 1709456 vanishing minors of order 4. For maxdet 
matrices of order n = 3 (mod 4), we sometimes find different numbers of 
vanishing minors of order 2. For example, if n = 11, we get 1391, 1389, and 
1401 vanishing minors for the three HT-equivalence classes. The right-hand- 
side of eqn. fflOl) . which by [9J Corollary 2] is a lower bound on the number 
of vanishing minors in this non-Hadamard case, gives 1362 (rounded up). 
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Finally, we briefly consider the frequencies (or multiplicities) with which 
the different values of | det(M)/2 n_1 | occur for minors of order m of a maxdet 
matrix of order n. In Table 5 we give the results of computations for m = 7, 
n — 15, which gives the typical behaviour that we have observed^ The sec- 
ond column gives the observed multiplicity of a minor with | det(M)|/2 n_1 
equal to the integer in the first column. The third column gives the multiplic- 
ities observed when taking a random sample of ( 1 ? 5 ) = 41409225 uniformly 
distributed {±l}-matrices of order m (we call this the random model). It 
is clear from the table that the actual distribution is nothing like the distri- 
bution for the random model. A \ 2 test gives an absurdly small probability 
< 1CP 10 that the two samples were drawn from the same distribution. Sim- 
ilar behaviour occurs for other values of m > 2. For example, when m = 2 
we find 5187 zero minors and 5838 nonzero minors, but for random matrices 
of order 2 we expect zero and nonzero values to occur with equal probability. 



|det(M)|/2 n " 1 | 


multiplicity 
of minors 


multiplicity in 
random model 


ratio 





12857784 


24030613 


0.54 


1 


8402100 


11140444 


0.75 


2 


10831128 


4662108 


2.32 


3 


3483909 


924336 


3.77 


4 


3935280 


504938 


7.79 


5 


622842 


76496 


8.14 


6 


927162 


55811 


16.61 


7 


129576 


7769 


16.68 


8 


201900 


6102 


33.09 


9 


17544 


608 


28.86 


total 


41409225 


41409225 





Table 5: Comparison of observed multiplicities of minors of 
order 7 in a maxdet matrix of order 15 with a random model 



Table 5 shows that the normalised minors are biased towards even val- 
ues. For the random model, this bias can be explained by reducing to the 
{0, 1} case and considering the evaluation of the determinant in Z/2Z. Then, 
for large n, we expect even values to occur about 71% of the timejfl This 

4 Data for other values of m and n may be found at [4]. 

5 The precise constant in the limit as n — > oo is 1 — rifc>i(l — 2~ fe ), see [51 UM 127] . 
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prediction is in agreement with the data for the random model (the third 
column). For the second column we find that even values occur about 69% 
of the time, which is close to the prediction for the random model. Thus, 
although the actual distribution differs considerably from that of the random 
model, the bias towards even values persists. In comparison with the random 
model, the extreme bias in favour of even values in the tail of the distribution 
compensates for a bias against zero minors. 

5 Algorithms for computing the set of minors 

Recall that Sk(A) is the set of k x k submatrices of an n x n {±l}-matrix A, 
where < k < n. In this section we describe algorithms for enumerating all 
the minors of A. Our application is to maxdet matrices, but the algorithms 
apply to all square {±l}-matrices A. 

In our enumeration we consider only the absolute values of minors, nor- 
malised by the factor 2™ _1 which always divides the determinant of an n x n 
{±l}-matrix, so the set of minors of A is defined to be the set 



possibilities to consider when finding the set S(A) for an n x n matrix A. 
In this section we consider four possible algorithms (of which we used two) 
for finding S(A). Their running times all involve the factor 4 n , so none of 
them is practical for n much larger than 20, but they differ significantly in 
the factor multiplying 4 n and in their space requirements. 

Each algorithm has two variants: the first just determines the set S(A) 
of minors; the second also counts the multiplicity of each minor, that is, how 
often a given value d G S(A) occurs. We describe the first variant of each 
algorithm, and briefly mention the changes required for the second variant. 

In the descriptions of the algorithms we explain how to compute the 
complete set S(A); it should be clear how to compute the subset of S(A) 
corresponding to the minors of a given order k. 





15 



Algorithm A 

Algorithm A simply considers, for each k in {1, 2, . . . , n}, the set of all k x k 
sub matrices of A, and evaluates the determinant of each such submatrix M by 
Gaussian elimination with partial pivoting, using floating-point arithmetic. 
The computed determinant is scaled by division by 2 n_1 and rounded to the 
nearest integer. 

Clearly there is a danger that rounding errors during the process of Gaus- 
sian elimination could lead to an incorrectly rounded integer result. However, 
our experiments, using IEEE standard 64-bit floating-point arithmetic PQ, 
showed that this is not a problem for the values of n that we considered 
(n < 25). Gaussian elimination with partial pivoting is numerically sta- 
ble [TBI |4"3"] . and the maximum scaled determinant of a 25 x 25 {±l}-matrix 
is 42 x 6 11 = 15237476352, meeting the Barba bound [31 [32], so only requires 
34 bits of precision, significantly fewer than the 53 bits provided by IEEE 
standard arithmetic. As a precaution, our implementation prints a warning 
and halts if the fractional part of a scaled determinant exceeds 1/8; this never 
occurred for n < 25. 

Gaussian elimination requires 0(k 3 ) arithmetic operations to evaluate the 
determinant of a k x k matrix. As is traditional in numerical analysis, we 
count multiplications/divisions but ignore additions/subtractions. With this 
convention, Gaussian elimination requires k 3 /3 + 0(k 2 ) operations. Thus the 
total cost is 



The storage requirements of Algorithm A are minimal, apart from the space 
required to store the results, that is the set S(A) of minors and (if required) 
their multiplicities. This is common to all the algorithms considered - they 
all need space to store their results. 

The set S(A) can be represented using one bit for each possible value 
| det |/2 n_1 . From the Hadamard bound, this requires at most 2 1 ~ n n n ^ 2 + 1 
bits. For example, if n = 20, it requires 19531251 bits (2.33 MB). For 
n = 24 it requires 519 MB, which is still feasible. For the variant that counts 
multiplicities, each bit needs to be replaced by an integer word (say 32 bits or 
4 bytes), so the storage required would be acceptable for n = 20 (75 MB) but 
excessive for n = 24 (16 GB), since the computers available to us typically 
have memories of 1 to 4 GB. 
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Fortunately, a much more economical representation of S(A) is usually 
possible, because not all minors in the range [0, L2 1_n n n//2 J] actually occur. 
The set S(A) is usually quite sparse, especially when the order n of A is 
divisible by 4. For example, with Hadamard matrices of order 16, we have 
#S(A) < 100 (see g}. Thus, instead of using 131073 bits to represent S(A), 
we can use a hash table with say 200 words [22] ■ With such an implementa- 
tion, the storage requirements are moderate for n < 25. 

Algorithm B 

Algorithm B is similar to Algorithm A, but uses a rank-1 updating formula 
to update the inverse and determinant of each k x k submatrix B of A if we 
already know the inverse and determinant of a submatrix that differs from 
B in only one row or column. The inverse updating formula 

(B + uv T )~ 1 = B- 1 - (1 + v T B- 1 u)- 1 B- 1 uv T B- 1 

is known at the Sherman-Morrison formula [3B] - the determinant updating 
formula 

det(B + uv T ) = det(£)(l + v T B- 1 u) 

seems to be known only as a "matrix determinant lemma" . 

Since the updating steps require ~ k 2 operations, the complexity is 

Wb ~ 4 n n 3 / 2 / y/n , 

so W a /Wb ~ n/24. 

We did not use Algorithm B because the constant factors involved make 
it slower than Algorithm A for n < 20, and because it is difficult to guarantee 
a correctly rounded integer result due to possible numerical instability. Also, 
even with exact arithmetic, we would have to use a different method whenever 
B is singular. 

Algorithm C 

Algorithm C uses integer arithmetic and evaluates each k x k determinant 
using the method (attributed to Laplace) of expansion by minors, see for 
example [301 Chapter 4]. To compute a k x k minor det(S), we need to 
know the (k — 1) x {k — 1) minors formed by deleting the first row and 
an arbitrary column of B. If these minors have been saved from a previous 
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computation, then the work involved in computing one minor | det(B) | is only 
k multiplications (by ±1) and k — 1 integer additions, plus any overheads 
involved in retrieving the previously stored values. If we assume, for purposes 
of comparison with Algorithms A-B, that the work involved amounts to k 
operations, then the total cost is 




giving Wa/Wc ~ n 2 /12. 

Unfortunately, this algorithm has a potentially large memory require- 
ment. If we compute the minors of order k in increasing order k — 1, 2, . . . , n, 
then to compute the minors of order k we need all minors of order 

k — 1. In the worst case, when k ~ n/2, the memory required to store the 
minors of order k — 1 and k is about 4 n+1 / (7™) words, which is too large to be 
practical for the values of n that we wish to consider. More memory-efficient 
implementations are possible, but complicated. For this reason we discarded 
Algorithm C and implemented a slightly slower, but much simpler algorithm, 
Algorithm D. 

Algorithm D 

The idea of Algorithm D is the same as that of Algorithm C. However, when 
computing the minors of order k of an nxn matrix A, the outer loop runs over 
all (™) combinations of k rows of A. Having selected these k rows, forming a 
kxn matrix B, we now compute all minors of order k of B. At the j-th step 
we compute all minors of order j in the last j rows of B. Thus the number 
of operations (counting as for Algorithm C) to compute the minors of order 
k of B is 




and the space requirement is at most 2(^ 2 j) ~ 2 n+ 2 f ^fwn words, much less 
than for Algorithm C. The overall operation count is 

k=i ^ 7 j=i vj/ 
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This is larger than the operation count for Algorithm C by a factor 
~ y/Tm/2, but smaller than the operation count for Algorithm A by a factor 
~n 3 / 2 /(6v^F). 

If a parallel implementation is desired, then it is easy to parallelise over 
the outer loop - different processors can work on different combinations of k 
rows of A in parallel. 

We ran both Algorithms A and D on small cases to check the correctness 
of our implementations. For the large cases we used mainly Algorithm D, 
which is much faster than Algorithm A for the most time-consuming cases 
{k ps n/2), as expected from the operation counts given above. 
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Appendix: Tables of minors for orders < 21 

In Tables 6-43, k — [n/A\ , where n is the order of the {±l}-matrix. The first 
column gives the order m of the minor, 1 < m < n. The second column gives 
the set of absolute values of the minors of order m, divided by the known 
factor 2 m_1 . In this column the notation "{a, b, . . .} x k a " is a shorthand for 
u {ak a , bk a , . . .}", etc. For n e {19, 21} we have abbreviated the entry in the 
second column by giving only the minimum and maximum rather than the 
complete set, using "(min, max)" instead of "{a, b, ...}". In such cases we 
write "(a, b) x k an instead of u {ak a ,bk a f. 

In the third column we give the scaled maximum determinant A(m) = 
D(m)/2 m ~ 1 (redundant, but included for easy comparison with the entries in 
the second column). The fourth column answers whether some minor meets 
the maximum possible determinant for its order (see Table 2), and the last 
column answers whether the full spectrum of possible values of minors (as 
given in Table 1) occurs. 

If there is more than one HT-equivalence class for an order n, the classes 
are listed in the same order as they are given in [32] . The information given 
in Tables 6-43 is sufficient to uniquely identify each HT-equivalence class. 

To avoid making this Appendix excessively long, we have omitted details 
such as the frequency of occurrence of each minor value. Further information 
is available on our website [1]. 
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Orders 1—6 



m 


{minors} 


A(m) 


max? 


full? 


1 


{1} 


1 


yes 


yes 


Table 6: n -- 


= 1, k = 


0, m,f = 


= 1 


m 


{minors} 


A(m) 


max? 


full? 


2 


{1} 


1 


yes 


no 


1 


{1} 


1 


yes 


yes 



Table 7: n = 2, k = 0, d = 1, m f = 1 



m 


{minors} 


A(m) 


max? 


full? 


3 




1 


yes 


no 


1-2 


full spectrum 


1 


yes 


yes 


Table 8: n = 3, k 


= 0,d = 


= 1, m f 


= 2 


m 


{minors} 


A(m) 


max? 


full? 


4 


{2} 


2 


yes 


no 


3 


{1} 


1 


yes 


no 


1-2 


full spectrum 


1 


yes 


yes 


Table 9: n = 4, k 


= l,d = 


= 1, m f 


= 2 


m 


{minors} 


A(m) 


max? 


full? 


5 


{3} 


3 


yes 


no 


4 


{1,2} 


2 


yes 


no 


1-3 


full spectrum 


1 


yes 


yes 


Table 10: n = 5,k 


= 1, d~- 


= 1, m f 


= 3 


m 


{minors} 


A(m) 


max? 


full? 


6 


{5} 


5 


yes 


no 


5 


{1-3} 


3 


yes 


no 


1-4 


full spectrum 


< 2 


yes 


yes 


Table 11: n = 6, fc 


= 1, d = 


= !, m f 


= 4 
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Orders 7-1 1(a) 



m 


{minors} 


Aim) 


max? 


full? 


7 


{9} 


9 


yes 


no 


1-6 


full spectrum 


< 5 


yes 


yes 


Table 12: n = 7, k 


= l,d = 


= 1, "V 


= 6 


m 


{minors} 


A(m) 


max? 


full? 


8 


{2} x 2 4 


2 x 2 4 


yes 


no 


7 


{1} x 2 3 


9 


no 


no 


6 


{0,1} x 2 2 


5 


no 


no 


5 


{0,1} x 2 1 


3 


no 


no 


1-4 


full spectrum 


< 2 


yes 


yes 



Table 13: n = 8, k = 2, d = 4, m f = 4 



m 


{minors} 


A(m) 


max? 


full? 


9 


{7} x 2 3 


7 x 2 3 


yes 


no 


8 


{2,3,4,6,8} x2 2 


8 x 2 2 


yes 


no 


7 


{0..4} x 2 1 


9 


no 


no 


1-6 


full spectrum 


< 5 


yes 


yes 



Table 14: n = 9, k = 2, d = 1, m f = 6 



m 


{minors} 


A(m) 


max? 


full? 


10 


{9} x 2 4 


9 x 2 4 


yes 


no 


9 


{3,6} x 2 3 


7 x 2 3 


no 


no 


8 


{0..5,8} x 2 2 


8 x 2 2 


yes 


no 


7 


{0..4} x 2 1 


9 


no 


no 


1-6 


full spectrum 


< 5 


yes 


yes 



Table 15: n = 10, k = 2, d = 2, m f = 6 



m 


{minors} 


A(m) 


max? 


full? 


11 


{20} x 2 4 


20 x 2 4 


yes 


no 


10 


{0,2..8,10,12,16} x 2 3 


18 x 2 3 


no 


no 


9 


{0,2..4,6,8..10,12..24, 
26.-33,36,40,48} 


56 


no 


no 


8 


{0..18,20,24} 


32 


no 


no 


1-7 


full spectrum 


< 9 


yes 


yes 



Table 16: n = 11(a), k = 2, d = 4, m f = 7 
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Orders ll(b)-12 



m 


{minors} 


A(m) 


max? 


full? 


11 


{20} x 2 4 


20 x 2 4 


yes 


no 


10 


{0, 1,4..6,8..11,14,16} x 2 3 


18 x 2 3 


no 


no 




{0,2..4,6,8..24,26..28, 


56 






9 


30..32,36,40,44,48} 


no 


no 


8 


{0..18,20,24} 


32 


no 


no 


1-7 


full spectrum 


< 9 


yes 


yes 



Table 17: n = 11(b), k = 2, d = 4, m f = 7 



m 


{minors} 


A(m) 


max? 


full? 


11 


{20} x 2 4 


20 x 2 4 


yes 


no 


10 


{4,6,8,12,16} x 2 3 


18 x 2 3 


no 


no 


9 


{0.-4,6,8,12} x 2 2 


14 x 2 2 


no 


no 


8 


{0.-8,12,16} x 2 1 


16 x 2 1 


yes 


no 


7 


{0-8} 


9 


no 


no 


1-6 


full spectrum 


< 5 


yes 


yes 


Table 18: n = 11(c), k 


= 2, d = 


3, rrif = 


6 



m 


{minors} 


A(m) 


max? 


full? 


12 


{2} x 3 6 


2x3 6 


yes 


no 


11 


{1} x 3 5 


5 x 2 6 


no 


no 


10 


{0,1} x 3 4 


3 2 x 2 4 


no 


no 


9 


{0,1} x 3 3 


7 x 2 3 


no 


no 


8 


{0-2} x 3 2 


2 5 


no 


no 


7 


{0-3} x 3 1 


3 2 


yes 


no 


1-6 


full spectrum 


< 5 


yes 


yes 



Table 19: n = 12, k = 3, d = 5, m f = 6 
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Orders 13 15 



m 


{minors} 


A(m) 


max? 


full? 


13 


{5} x 3 6 


5 x 3 6 


yes 


no 


12 


{2,3} x 3 5 


6 x 3 5 


no 


no 


11 


{0..3} x 3 4 


5 x 2 6 


no 


no 


10 


{0-4} x 3 3 


3 2 x 2 4 


no 


no 


9 


{0..5} x 3 2 


7 x 2 3 


no 


no 


8 


{0-6} x 3 1 


2 5 


no 


no 


1-7 


full spectrum 


< 9 


yes 


yes 



Table 20: n = 13, k = 3, d = 6, m f = 7 



m 


{minors} 


A(m) 


max? 


full? 


14 


{13} x 3 6 


13 x 3 6 


yes 


no 


13 


{4,6,7,9} x 3 5 


15 x 3 5 


no 


no 


12 


{0..7,9,10} x 3 4 


18 x 3 4 


no 


no 


11 


{0..9, 11} x 3 3 


5 x 2 6 


no 


no 


10 


{0..13} x 3 2 


3 2 x 2 4 


no 


no 


9 


{0-15} x 3 1 


7 x 2 3 


no 


no 


8 


{0..18,20} 


2 5 


no 


no 


1-7 


full spectrum 


< 9 


yes 


yes 


Table 21: n = 14, k 


= 3, d = 


7, m f = 


7 



m 

14 
13 
12 
11 
10 

9 

8 

1-7 



{minors} 
{35} x 3 6 
{7, 8, 12, 14, 17, 18, 21, 23, 27} x 3 5 
{0.-21,23,24,26,27} x 3 4 
{0.. 22, 24.. 27} x 3 3 
{0..29,31,35} x 3 2 
{0..36,39,40} x 3 1 
{0..36, 38..40,42,44, 45} 
{0-18,20,24} 
full spectrum 



A(m) 


max? 


full? 


35 x 3 6 


yes 


no 


39 x 3 5 


no 


no 


45 x 3 4 


no 


no 


54 x 3 3 


no 


no 


5 x 2 6 


no 


no 


3 2 x 2 4 


no 


no 


56 


no 


no 


32 


no 


no 


< 9 


yes 


yes 



Table 22: n = 15, k = 3, d = 8, m f = 7 
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Order 16 



m 


{minors} 




max: 


lull : 


16 




{2} x 4 8 


2x4 


yes 


no 






{1} x 4 7 


OO X O 


no 


no 


14 




{0,1} x 4 6 


13 x 3 6 


no 


no 


13 




{0,1} x 4 5 


r. v 06 
X O 


no 


no 


12 




{0..2} x 4 4 


Z A O 


no 


no 


11 




{0..3} x 4 3 


X Z 


no 


no 


10 




{0..2} x 2 5 


Q v 9 4 

y x z 


no 


no 


9 




{0..2} x 2 4 


7 v 9 3 

1 A Z 


no 


no 


8 




{0..4} x 2 3 


^9 
oz 


yes 


no 


7 




{0..2} x 2 2 


Q 


no 


no 


6 




{0..2} x 2 1 




no 


no 


1-5 


full spectrum 


< 3 


yes 


yes 


Table 23: n = 16(a), 


k = 4, d 


= 8, rrif = 5 


m 


{minors} 


A(m) 


max? 


full? 


11-16 


as for 16(a) 






no 


10 




{0..5} x 4 2 


9 x 4 2 


no 


no 


9 




{0..4} x 2 3 


7x 2 3 


no 


no 


8 




{0..6,8} x 2 2 


32 


yes 


no 


7 




{0..4} x 2 1 


9 


no 


no 


1-6 




full spectrum 


< 5 


yes 


yes 



Table 24: n = 16(b), k -■ 


= 4, d = 8 


,, m/ = 


6 


m 


{minors} 


A(m) 


max? 


full? 


10-16 


as for 16(b) 






no 


9 


{0..9} x 4 1 


14 x 4 1 


no 


no 


8 


{0.-10,12,16} x 2 1 


16 x 2 1 


yes 


no 


1-7 


full spectrum 


< 9 


yes 


yes 


Table 25: n = 16(c), k -- 


= 4, d = 8 


, m f = 


7 



m 


{minors} 


A(m) 


max? 


full? 


9-16 


as for 16(c) 






no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 26: n = 16(d), k = 4, d = 8, m/ = 8 
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Order 17 



m 


{minors} 


A(m) 


max? 


full? 


17 


{5} x 4 s 


5 x4 s 


yes 


no 


16 


{2,3,8}x4 7 


8 x 4 7 


yes 


no 


15 


{0..4} x 4 6 


35 x 3 6 


no 


no 


14 


{0-4} x 4 5 


13 x 3 6 


no 


no 


13 


{0..7} x 4 4 


5 x 3 6 


no 


no 


12 


{0-9} x 4 3 


2 x 3 6 


no 


no 


11 


{0..13,15} x 4 2 


20 x 4 2 


no 


no 


10 


{0..21,24,27} x 4 1 


36 x 4 1 


no 


no 


9 


{0.-40,42,44,45,48} 


56 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 27: n = 17(a), k = 4, d = 1, m f = 8 



m 


{minors} 


A(m) 


max? 


full? 


10-17 


as for 17(a) 






no 


9 


{0-22,24} x 2 1 


28 x 2 1 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 28: n = 17(b), k = 4, d = 1, m f = 8 



m 


{minors} 


A(m) 


max? 


full? 


11-17 


as for 17(a) 






no 


10 


{0-10,12} x 2 3 


18 x 2 3 


no 


no 


9 


{0-12} x 2 2 


14 x 2 2 


no 


no 


8 


{0-10,12,16} x 2 1 


16 x 2 1 


yes 


no 


1-7 


full spectrum 


< 9 


yes 


yes 



Table 29: n = 17(c), k = 4, d = 1, m f = 7 
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Order 18 



m 


{minors} 


A(m) 


max? 


full? 


18 


{17} x 4 8 


17 x 4 8 


yes 


no 


17 


{6,7, 10,11} x 4 7 


20 x 4 7 


no 


no 


16 


{0..8, 10, 11, 13} x 4 6 


32 x 4 6 


no 


no 


15 


{0..16} x 4 5 


35 x 3 6 


no 


no 


14 


{0-20} x 4 4 


13 x 3 6 


no 


no 


13 


{0..24,26,28} x 4 3 


5 x 3 6 


no 


no 


12 


{0.-38,40,41,44,52} x 4 2 


2 x 3 6 


no 


no 


11 


{0..62,64,68,80} x 4 1 


80 x 4 1 


yes 


no 


10 


{0..94,96..98, 100. .102, 


144 






104,108,112,128} 


no 


no 


9 


{0.-40,42,44,45,48} 


56 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 30: n = 18(a), k = 4, d = 7, m f = 8 



m 


{minors} 


A(m) 


max? 


full? 


14-18 


as for 18(a) 






no 


13 


{0.-24,26,28,32} x 4 3 


5 x 3 6 


no 


no 


12 


{0..38,40,41,44, 64} x 4 2 


2 x 3 6 


no 


no 


11 


{0-62,64,68,80} x 4 1 


80 x 4 1 


yes 


no 


10 


{0-52,54,56,64} x 2 1 


72 x 2 1 


no 


no 


9 


{0-40,42,44,45,48} 


56 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 31: n = 18(b), k = 4, d = 7, m f = 8 



m 


{minors} 


A(m) 


max? 


full? 


14-18 


as for 18(a) 






no 


13 


{0-24,32} x 4 3 


5 x 3 6 


no 


no 


12 


{0-37,40,64} x 4 2 


2 x 3 6 


no 


no 


11 


{0-32} x 2 3 


40 x 2 3 


no 


no 


10 


{0-28,32} x 2 2 


36 x 2 2 


no 


no 


9 


{0-22,24} x 2 1 


28 x 2 1 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 32: n = 18(c), k = 4, d = 10, m f = 8 
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Order 19 



m 


(min, max) minor 


A(m) 


max? 


full? 


19 


(833,833) x 4 y 


833 x 4 6 


yes 


no 


18 


(140, 784) x 4 5 


1088 x 4 5 


no 


no 


17 


(0, 672) x 4 4 


1280 x 4 4 


no 


no 


16 


(0,676) x 4 3 


2048 x 4 3 


no 


no 


15 


(0, 1050) x 4 2 


35 x 3 6 


no 


no 


14 


(0, 1470) x 4 1 


13 x 3 6 


no 


no 


13 


(0, 1904) 


3645 


no 


no 


12 


(0, 756) 


1458 


no 


no 


11 


(0,312) 


320 


no 


no 


10 


(0, 128) 


144 


no 


no 


1-9 


full spectrum 


< 56 


yes 


yes 


Table 33: n = 19(a), k = 4, d = 10, m f = 


o 


m 


(min, max) minor 


A(m) 


max? 


full? 


19 


(833,833) x 4 6 


833 x 4 6 


yes 


no 


18 


(168,616) x 4 5 


1088 x 4 5 


no 


no 


17 


(0, 672) x 4 4 


1280 x 4 4 


no 


no 


16 


(0, 740) x 4 3 


2048 x 4 3 


no 


no 


15 


(0, 1024) x 4 2 


35 x 3 6 


no 


no 


14 


(0, 1536) x 4 1 


13 x 3 6 


no 


no 


13 


(0, 2048) 


3645 


no 


no 


12 


(0, 1024) 


1458 


no 


no 


11 


(0, 288) 


320 


no 


no 


10 


(0, 128) 


144 


no 


no 


1-9 


full spectrum 


< 56 


yes 


yes 



Table 34: n = 19(b), k = 4, d = 10, m f = 9 



m 


(min, max) minor 


A(m) 


max? 


full? 


14-19 


as for 19(b) 








13 


(0,2560) 


3645 


no 


no 


1-12 


as for 19(b) 









Table 35: n = 19(c), k = 4, d = 10, m f = 9 
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Order 20 



m 


{minors} 


A(m) 


max? 


full? 


20 


{2} x 5 1U 


2 x 5 10 


yes 


no 


19 


{1} x 5 9 


833 x 4 6 


no 


no 


18 


{0,1} x 5 8 


17 x 4 8 


no 


no 


17 


{0,1} x 5 7 


5 x 4 8 


no 


no 


16 


{0..2} x 5 6 


2 x 4 8 


no 


no 


15 


{0-3} x 5 5 


35 x 3 6 


no 


no 


14 


{0-5} x 5 4 


13 x 3 6 


no 


no 


13 


{0-9} x 5 3 


5 x 3 6 


no 


no 


12 


{0.-18,20,24,32} x 5 2 


2 x 3 6 


no 


no 


11 


{0.-40,42,44,48} x 5 1 


64 x 5 1 


no 


no 


10 


{0..92, 95, 96, 100, 104, 108, 112, 125, 144} 


144 


yes 


no 


9 


{0.-40,42,44,48} 


56 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 



Table 36: n = 20(a), k = 5, d = 10, m f = 8 



m 


{minors} 


A(m) 


max? 


full? 


12-20 


as for 20(a) 






no 


11 


{0.-40,42,44,45,48} x 5 1 


64 x 5 1 


no 


no 


10 


{0..90,92,93,95..102,104,108, 


144 






112,117,120,125,128,144} 


yes 


no 


9 


{0-40,42,44,45,48} 


56 


no 


no 


1-8 


full spectrum 


< 32 


yes 


yes 




Table 37: n = 20(b), k = 5, d 


= 10, rrif 


= 8 





m 


{minors} 


A(m) 


max? 


full? 


11-20 

10 

1-9 


as for 20(b) 
{0..88, 90, 92, 93, 96, 99, 100, 102, 
104, 108, 112, 120, 125, 128, 144} 

as for 20(b) 


144 


yes 


no 
no 



Table 38: n = 20(c), k = 5, d = 10, m f = 8 



32 



Order 21(a)-(b) 



m 


(min, max) minor 


A(m) 


max? 


full? 


21 




(29,29) x 5 9 


29 x 5 9 


yes 


no 


20 




(10,30) x 5 8 


50 x 5 s 


no 


no 


19 




(0,35) x 5 7 


833 x 4 6 


no 


no 


18 




(0,40) x 5 6 


17 x 4 s 


no 


no 


17 




(0,45) x 5 5 


5 x 4 8 


no 


no 


16 




(0,65) x 5 4 


2 x 4 8 


no 


no 


15 




(0, 100) x 5 3 


35 x 3 6 


no 


no 


14 




(0, 240) x 5 2 


13 x 3 6 


no 


no 


13 




(0,416) x 5 1 


5x3 6 


no 


no 


12 




(0,800) 


1458 


no 


no 


11 




(0,320) 


320 


yes 


no 


10 




(0, 144) 


144 


yes 


no 


1-9 




full spectrum 


< 56 


yes 


yes 


Table 39: n = 21(a), k 




- 5, d = 10, rrif = 


9 


m 


(min, max) minor 




A(m) 


max? 


full? 


15-21 


as for 21(a) 






no 


14 




(0,216) x 5 2 


13 x 3 6 


no 


no 


13 




(0,400) x 5 1 


5 x 3 6 


no 


no 


12 




(0, 800) 


1458 


no 


no 


11 




(0,320) 


320 


yes 


no 


1-10 




full spectrum 


< 144 


yes 


yes 



Table 40: n = 21(b), k = 5, d = 10, m f = 10 
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Order 21(c)-(e) 



m 


(min, max) minor 


A(m) 


max? 


full? 


15-21 


as for 21(a) 






no 


14 


(0,200) x 5 2 


13 x 3 6 


no 


no 


13 


(0,400) x 5 1 


5 x 3 6 


no 


no 


12 


(0,800) 


1458 


no 


no 


11 


(0, 304) 


320 


no 


no 


10 


(0, 144) 


144 


yes 


no 


1-9 


full spectrum 


< 56 


yes 


yes 


Table 41: n = 21(c), k = 


5, d — 11, rrif = 


Q 


m 


(min, max) minor 


A(m) 


max? 


full? 


15-21 


as for 21(a) 






no 


14 


(0,240) x 5 2 


13 x 3 6 


no 


no 


13 


(0,416) x 5 1 


5 x 3 6 


no 


no 


12 


(0, 800) 


1458 


no 


no 


11 


(0,320) 


320 


yes 


no 


1-10 


full spectrum 


< 144 


yes 


yes 


Table 42: n = 21(d), k = 


5, d — 10, rrif = 


10 


m 


(min, max) minor 


A(m) 


max? 


full? 


15-21 


as for 21(a) 






no 


14 


(0,212) x 5 2 


13 x 3 6 


no 


no 


13 


(0,368) x 5 1 


5 x 3 6 


no 


no 


12 


(0, 800) 


1458 


no 


no 


11 


(0,288) 


320 


no 


no 


10 


(0, 144) 


144 


yes 


no 


1-9 


full spectrum 


< 56 


yes 


yes 



Table 43: n = 21(e), k = 5, d = 11, m f = 9 
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